The effects of non-abelian statistics on two-terminal shot noise in a quantum Hall 

liquid in the Pfaffian state 
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We study non-equilibrium noise in the tunnelling current between the edges of a quantum Hall 
liquid in the Pfaffian state, which is a strong candidate for the plateau at v = 5/2. To first non- 
vanishing order in perturbation theory (in the tunneling amplitude) we find that one can extract 
the value of the fractional charge of the tunnelling quasiparticles. We note however that no direct 
information about non-abelian statistics can be retrieved at this level. If we go to higher-order in 
the perturbative calculation of the non-equilibrium shot noise, we find effects due to non-Abelian 
statistics. They are subtle, but eventually may have an experimental signature on the frequency 
dependent shot noise. We suggest how multi-terminal noise measurements might yield a more 
dramatic signature of non-Abelian statistics and develop some of the relevant formalism. 

PACS numbers: 73.23.-b, 71.10.Pm, 73.43.Jn 



I. INTRODUCTION 

Recently, there has been increased interest in the pos- 
sibility that non-Abelian braiding statistics might occur 
in the fractional quantum Hall regime, particularly at 
v = 5/2 and, possibly, other fractions in the first excited 
Landau level. One possible state describing the v = 5/2 
plateau is the Pfaffian stateiSi 3 ^, which combines as- 
pects of the fractional quantum Hall effect and BCS pair- 
ing. Its quasiparticle excitations obey non-Abelian braid- 
ing statisticsi*2iii4. Thus far, the strongest evidence in 
favor of this state comes from numerical studies of elec- 
trons in the first excited Landau levels. Given the special 
topological properties of this state, it would be of great 
importance to have a clean experimental test to identify 
it. Such tests, measuring directly the bra idin g of quasi- 
particles. have been proposed in refs. RHlU lTM Here, we 
see if more indirect, but simpler experiments can reveal 
some of the same properties of this state. 

Current noise measurements are an interesting possi- 
bility. Since the current is delivered in fractional charge 
e* packets, the zero-frequency limit of the current noise 
evinces shot noise. The ratio between the shot noise 
and the current is simply e*, as has been observed 
experimentally^. It is believedii*^ that such experi- 
ments could also be used in various setups to also "mea- 
sure" the statistics of quasiparticles, though actual ex- 
periments addressing this issue are much fewer—. In 
photon-counting experiments, "bunching" is observed, 
while the analogous noise measurments in a Fermi liquid 
show "anti-bunching" , as a result of Fermi statistics. In 
a state in which the quasiparticles are anyons, one might 
expect something intermediate between these two lim- 
its, from which the statistics might be extracted. In the 
Pfaffian state, something even more interesting might be 
observed since the quasiparticles have non-Abelian statis- 
tics. 

The difference between a non-Abelian state, such as 



the Pfaffian, and an Abelian state, such as the Laugh- 
lin state, is most easily seen when two quasiparticles are 
brought together, or 'fused'. In an Abelian state, there is 
a unique quasiparticle type which is obtained (up to addi- 
tional bosonic excitations). For instance, when two e/3 
Laughlin quasiparticles are brought together, a charge 
2e/3 quasiparticle is obtained. When a charge e/3 quasi- 
particle is taken around a 2e/3 quasiparticle, the wave- 
function acquires a factor e 47 ™/ 3 . One the other hand, 
when two charge e/4 quasiparticles in the Pfaffian state 
are fused together, there are two possible outcomes, both 
with charge e/2. They differ by the presence or absence 
of a neutral fermion. Thus, when a charge e/4 quasi- 
particle is taken around two fused e/4 quasiparticles, the 
factor which is acquired is either i or —i, depending on 
how the particles fuse. (When they fuse to form a linear 
combination of the two possible outcomes, taking an e/4 
quasiparticle around them would transform them into a 
different linear combination.) It is this multiplicity of 
outcomes of fusion which we would like to access through 
consideration of the noise. 

In this paper we analyze the non-equilibrium noise in 
the tunnelling current between two edges of a FQH liq- 
uid in a Pfaffian state and we find that, indeed, infor- 
mation about the fractional charge and statistics of the 
FQH Pfaffian quasiparticles can be extracted from the 
two-terminal noise. In particular, we find that up to first 
order in perturbation theory, the shot noise contains suf- 
ficient information to allow for the measurement of the 
fractional charge of the tunnelling quasiparticles. The 
charge also shows up in the Josephson frequency. When 
we also compute higher-order corrections to the noise, 
non-Abelian effects come into play. One of the effects 
of statistics in the Pfaffian state is that the peak in the 
noise at the Josephson frequency is shifted toward lower 
frequencies and enhanced, while in the Laughlin state it 
is flattened. 

This is less dramatic than the underlying physics be- 
cause once the quasiparticle four-point function correla- 
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tion function is Fourier transformed and inserted into the 
expression for the noise, its structure is masked. Nev- 
ertheless, a two-terminal noise measurement is a good 
starting point for discussing more complicated setups. 
We discuss other possible experiments, with more termi- 
nals, and assess their usefulness for observing the effects 
of non-Abelian statistics. 

In Section ^ we lay out the formalism necessary to 
calculate the tunnelling current and current fluctuations 
and note the expected non-linear I — V characteristic 
for tunneling into the Pfaffian state. In Section III we 
analyze the first non-zero order in perturbation theory 
for both the tunnelling current and shot noise. In Section 
IV we go to the next order in perturbation theory and to 
finite temperature. In Section V we present and discuss 
our results. We conclude in Section VI. 



II. FORMALISM 

The Pfaffian statei*2i2i4 describes electrons in a half- 
filled Landau level (with straightforward generaliza- 
tion to other even-denominator filling fractions or odd- 
denominator filling- fractions for systems of bosons). 
While there does not appear to be a FQH plateau at 
v = 1/2, there is one at v = 5/2 = 2 + 5. We ignore 
the filled lowest Landau level (of both spins), and fo- 
cus on the half-filled first excited Landau level, which we 
suppose is described by the Pfaffian state. The gapless 
chiral theory describing excitations at the edge of a Pfaf- 
fian state is a 1 + 1 — D conformal field theory (CFT) 9 . 
The edge theory has a charge sector - which is a free 
boson - and a neutral sector which is a c = 1/2 Majo- 
rana fermion. The latter is the chiral part of the critical 
theory of the 2D Ising model. In addition to the iden- 
tity operator, the latter has a spin or twist field a and a 
Majorana fermion operator -0. 

Considering the fact that a physical operator must 
have a single-valued correlation function with the elec- 
tron operator, the operators which can be identified with 
quasiparticle operators are 

$ 1/2 (x) = eV^*); $i(x) = e^ Hx) ; the neutral 

fermion ip; and \t e (x) = rpe t '^^ x \ Here, g = 1/2. 
Here $ e is the electron annihilation operator, while 
$1/4, and $1/2 annihilate quasiparticlcs with fractional 
charge ge/2 — e/4, and ge = e/2 respectively. In 
a tunnelling process between two edges of a quantum 
Hall liquid, the most relevant (in the RG sense) pro- 
cess will be the backscattering of the e/4 charged $1/4 
quasiparticles. We need to note also that the fields 
</>, ip, and a also need to be characterized by a chi- 
ral index R/L, denoting which of the right/left moving 
branches of the theory the fields belong to, such that 
*i/*,r/l(x) = e ±l ^*'^l 2 a R/L {x), and 9% /L (x) = 

We can write an edge theory for the charged cj> sector 
of the theory, which is analogous to the Laughlin state 



with g = 1/2. The corresponding Lagrangian density for 
the right and left moving modes isi^*i&: 

C R/l = ^dx<pR/ L (±dt + vd x )<p R/L , (1) 

where v is the velocity of the edge excitations, which 
we will set to 1. The fields 4>r/l satisfy the equal time 
commutation relations 

[<Pr/l (t,x), <j>B,/L (t,y)}= ±iTrsgn(x -y). (2) 

We can also writ o 15 ' 16 the total Lagrangian in terms of 

4> = 4>R + <t>L, 

Co(x,t) = i{[d t #M)] 2 - [dAM] 2 }- (3) 

where <fr is satisfying [<j>(t, x), d t <j>(t, y)] = 4:iri5(x — y). 
The Lagrangian for the neutral sector of the theory is 

£r/l = ^ ^R/h{±dt + vd x )ip R/L , (4) 

The Ising spin fields <j r /l are twist fields for the Majo- 
rana fermion ip R /L : whenever a Majorana fermion goes 
around its twist field, it changes sign. The Ising spin 
are most easily described in a conformal field theory pic- 
ture. The chiral ctr/i, fields have conformal weight 1/16, 
while the combination a = cj r <jl, which we will use later 
on, has dimension 1/8. Also, information about the two 
point correlations and four point correlation functions of 
the a fields have been derived in the CFT context. While 
the two-point functions are simple power laws, the four 
point functions have a more complicated structure, which 
incorporates the effects of non-Abelian statistics. This is 
described in detail in Appendix B. 

We are going to study a FQH setup similar to the one 
depicted in FigQ] In the presence of an applied voltage 
V, a current Iq is injected in the sample through lead 
A. In the absence of inter-edge tunnelling, this current 
will be picked up in lead C. If a gate voltage is applied 
across the sample, such that the sample is constricted, 
as indicated in the figure, quasiparticlcs from one edge 
can tunnel to the other edge, thereby giving rise to a 
tunnelling current It- This current can be measured in 
lead B. While the current in lead A will not be changed, 
the current picked up in lead C will be reduced to Iq — I t . 

We can write down the tunneling operator between 
the two edges of the quantum Hall fluid noting that a 
tunnelling process annihilates a right/left mover quasi- 
particle and creates a left/right moving quasiparticle on 
the opposite moving branch (see Fig^). 

#tu„ = r$ f lML (a; = 0)$ 1/4 , fi (.T = 0) + h.c, (5) 

which generates a change in the Lagrangian density: 

5L(t,x) = -Ta R (t,0)a L (t,0)e i ^ t ' ^/ 2 + h.c. (6) 

From the scaling dimensions of the fields in J5J, we de- 
duce that the tunneling operator has dimension 3/4. 
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In the next section we will find the expectation value of 
the tunnelling current operator and of the shot noise up 
to the first non-zero order in perturbation theory in the 
tunnelling amplitude T. 



FIG. 1: The two-terminal noise setup in a quantum Hall bar. 
Current is injected at A and measured at B. The voltage 
drop between these two terminals is also measured. 



III. NON-EQUILIBRIUM SHOT NOISE AND 
TUNNELLING CURRENT TO FIRST NON-ZERO 
ORDER IN PERTURBATION THEORY 



Hence, we expect the non-linear I-V characteristic for 
weak inter-edge tunneling at T = to bei^ 

It ~ V- 1 ' 2 . (7) 

For T > V, the tunneling conductance varies as 

Gt ~ T- 3 / 2 . (8) 

Note that the operator which tunnels a charge- 1/2 
quasiparticle is also relevant, but less so than JSJ. 

#tun = r'$ f 1/2 L {X = 0)$ 1/2 , R ( X = °) + h-C. 

= _r e ^ ( '*< 0) v^ + h.c. (9) 

A voltage drop between the edges of the quantum Hall 
liquid can be introduced 15 i 16 by letting T — > Te~ luJot , 
where ujq = luj = geV/2h. The tunnelling current op- 
erator is I t = ^[N R ,H] = -^[Nl.H] where N R/L are 
the total charge operators on the R/L edges. Using the 
commutation relations between the charge and the quasi- 
particle operators^, we find 

I t (t) ^ie*Tal(t,0)a R (t,0)e^ t ^^/ 2 e~ ZbJot + h.c, (10) 

where e* = ge/2 — e/4. The fluctuations (shot noise) in 
the tunnelling current can be written as 

S(u) = J Stye™, (11) 

where 

S(t) = l{I t (t),I t (0)}. (12) 
I 



We note that our problem is non-equilibrium, so that 
we should use the Keldysh non-equilibrium formalism 20 . 
However, for the zeroth and first order perturbation the- 
ory, there is no difference between the equilibrium and 
non-equilibrium formalisms^, so in this section we will 
make use of the equilibrium formalism to study the ex- 
pectation values of the current and shot noise to first 
non-zero order in perturbation theory. Thus 

(It) = \ J Ite lS , (13) 

where S is the total action S = So + SS = J dt£o + 
j dtSC, and Z is the partition function. We can expand 
the exponential to first order in the tunnelling to obtain 



(I t )=i dt'{I t (t)SC(t')) , (14) 



where ()o denotes taking expectation values with respect 
to the non-perturbed action Sq. 

Since the expectation value of the current It is time 
independent, we can set t — 0. Also, one cannot dis- 
tinguish between cr R ^ L and <j r /li so we wm write our 
result in terms of the expectation value of the operator 
(7 = <jr(Jl. Moreover, since x — for all operators, we 
will drop the spatial index, and we will only write down 
the time index. Also, all expectation value symbols on 
the right hand side refer now to the unperturbed action. 
The expectation value for the current becomes 



(I t ) = e*\T\ 2 j A'(a(0)a(t'))[e i " ot '(e^ (0)/2 e" l ^ (t ' )/2 ) - e ~**>* (e-WftW/teWSH*)/*)]. (15) 

I 

We note that ( e ±Vs<«t)/2 e =FV5<«0)/2) = e g{4>(tMo))/± that (a(z)a(z')) = \z - z'\- & ° , where 5 a = 1/4, and 
andi^ that (4>(t)<p(0)} = — 21n(e + it), where e is a short z = r + ix. This yields by analytical continuation 
time/high energy cutoff. Also, we know from CFT 21 , t — > it + e, (a(t)a(0)) = (e + it)~ s " . 
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Putting all the factors together and performing the 
integrals over time one obtains 

(It) = ^e*|lf^-\ (16) 

where is the Euler Gamma function, S = 8 a + 5^ — 
S a + g/2 = 1/2, e* = ge/2 = e/4, and w = geV/2h 
is the Josephson frequency corresponding to the applied 
voltage V. 

Note that the total current injected in the sample 
through lead A due to the applied voltage V is Iq — 
(g / 2) (e 2 / h)V , which is splitting between the leads B and 
C, (It) being the tunnelling current picked up in lead B. 
Obviously, in the absence of tunnelling (T — 0) the tun- 
nelling current is zero and the currents in leads A and C 
are equal, and equal to Iq- 

Similarly one can calculate the shot noise in the tun- 
nelling current. Up to the first non-zero order in pertur- 
bation theory this is: 

(S(t)) =~<I(i)I(0)+/(0)J(t)) 

= e* 2 |r| 2 cos(w <)(CT(t)CT(0))e 9 ^( t ^(°^/ 4 + (t -» -t) 
= e-|r| 2 cos(^)[ I -^ F + I -^ F ]. (17) 

We note that in the limit of e 0, we can rewrite the 
expression for (S(t)) as follows. 

(S(t)) = e* 2 \T\ 2 cos(wo0 |4-[e l7r5s9 " (t)/2 + e -" 5s9n{t)/2 ] 

= 2e* 2 |r| 2 cos(woi) A cos(tt(5/2). (18) 
r I 

So (S(t)), even to first non-zero order in perturbation 
theory incorporates in it information about both the frac- 
tional charge and the fractional statistics of the quasi- 
particles. For example, e* = ge/2 is a measure of the 
fractional charge of the tunnelling quasiparticles. The 
factor cos(7r5/2), where 8 = g/2 + S a is a measure of the 
statistical angle, which incorporates combined informa- 
tion about the charge mode through g/2, as well as about 
the neutral Ising mode through S a . Also (S(t)) decays 
as a power law in time, with a coefficient S. Moreover, 
it oscillates with a periodicity ujq — geV/2h, the Joseph- 
son frequency, which incorporates information about the 
fractional charge ge/2. However, no direct information 
is retrieved about the non-abelian a Ising mode, since its 
non-abelian properties only appear for fourth or higher 
order correlation functions of the a operators^!. 

Note that for the simple case of tunnelling between the 
edge states of a Laughlin quantum Hall liquid with filling 
fraction our results change such that 6 — > Si = 2v, and 
e* —>■ e'\ = ve. 

We can also study the frequency dependent noise 
(S(oS)) = J t e lwt (S(t)). By taking a non-zero cutoff e, 



and performing the integral over t exactly, one finds that 

(S(u)) = ^e^imiu;-^] 5 - 1 + l^ + col 5 - 1 ] 

= ^-I t [\l - uj/cvof- 1 + \1 + Lu/cjof- 1 }- (19) 

We note that the ratio (S(ui)) / (It) is a universal function, 
independent of the strength of the tunnelling T. More- 
over, at zero frequency (S(uj))/(It) — e*, from where a 
precise measurements of the fractional charge can be per- 
formed, and indeed has been doneifi for the simpler case 
of a Laughlin state with e* = e/3. For our situation the 
fractional charge e* = ge/2 — e/4 of the tunnelling quasi- 
particles can be extracted. The last feature we note is the 
presence of the singularity in frequency at the Josephson 
frequency lu = ujq. 

The quantity (S(lu)) seems to yield less information 
than (S(t)), however, it allows for a clean measurement 
of various quantities such as the fractional charge and the 
Josephson frequency, independently of the magnitude of 
the tunnelling coefficient T or of other experimental fac- 
tors. We note that at this order in perturbation theory, 
it is impossible to extract direct information about the 
statistics of the non-abelian Ising mode. Since we are in- 
terested in the non-abelian characteristics of the Pfaffian 
state, we will go to the next order in perturbation the- 
ory, from where more information about the non-abelian 
character of the quasiparticles can be extracted. 

We should note that all our calculations till now are 
performed at zero temperature. The generalization to 
finite temperature is straightforward, but at this order 
in perturbation theory the finite temperature will not 
provide any novel information, so we will not study this 
situation until later on in the following section. 



IV. HIGHER ORDER PERTURB ATIVE 
CORRECTIONS TO SHOT NOISE 

In the previous section we computed the tunnelling 
current and the shot noise up to the first non-zero or- 
der in perturbation theory. In this section we will study 
the next order correction. Before starting the calcula- 
tion, a few comments are in order. First, we need to note 
that to this perturbative order, one also needs to sub- 
tract the noncorrelated piece (I t (t)) (I t {0)) = (It) 2 from 
the expectation value of the shot noise (S(t)). Subtract- 
ing this contribution has no effect to the order |r| 2 , as 
(It) 2 is of order |r| 4 , but needs to be considered when 
we compute the shot noise to the order |T| 4 . We expect 
that (S(t)) — (I t ) 2 goes to zero when t is large, as the 
current should be uncorrelated with itself for large time 
separations. Moreover, the tunnelling current (I t ) is in- 
dependent of time, so the correction to (S(t)) will be a 
mere constant. In principle this correction can be de- 
termined easily from our calculations, but since it may 
slightly depend on various cutoffs, etc., we will account 
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for it by automatically renormalizing (S(t)) such that it 
goes to zero for large t. 

The second thing we need to note is that we are usu- 
ally interested not in the actual (S(u>)), but in the ratio 
(S(co))/(I t ). If we compute (S(ui)) up to |T| 4 , in order 
to compute the ratio accurately to order |F| , then we 
should also carefully compute (It) to order |T| 4 . How- 
ever, since (It) is a constant with respect to the frequency 
uj, this will only renormalize the ratio between the first 
and the second order corrections in (S (uj)) / (I t ) , and thus 
could be easily accounted by a slight renormalization in 
the tunnelling coefficient T. For these reasons, in this 
paper we will neglect these higher order corrections to 
(h). 

However, we should note that (It) depends on the ap- 
plied voltage, so that if we want to consider the depen- 
dence of (S (cu)) / (It) at a fixed frequency u> on the ap- 
plied voltage, one needs to take also into account the 
corrections of order |T| 4 to (I t ). It actually may be a 
good exercise with interesting consequences on its own, 
to compute these higher order corrections in (I t ) and see 
what, if any, are the consequences of the non-abelian and 
fractional statistics for the tunnelling current. 

In order to compute the next order in perturbation 
theory for (S(t)), it does no longer suffice to resort to 
an equilibrium approach, we need to use a Schwinger- 
Keldysh non-equilibrium approach 20 . We thus double 
the integration contour such that our time integrals ex- 
tend from — oo to oo (the r\ = + part of the contour), and 
back from oo to — oo (the r\ = + part of the contour) (see 
Fig|2l. We thus introduce a new index ?y to our opera- 
tors which describes the branch of the contour on which 
a specific operator sits. In terms of the new operators, 




— o=^ + axj 



FIG. 2: The two branches 77 = ± of the Keldysh contour 
we can write 

s(t) = l[i+(t)i-(o) + i+(o)r(t)} (20) 

and 

w*)) = \ E^w-'W^'^o, (21) 

where again we suppressed the spatial indices and the 
()o denotes expectation value with respect to the free 
action J K dt J d 3 xCo. The symbol Tk denotes time or- 
dering along the Keldysh contour described above and 
depicted in Fig. and /„ denotes a time integral 
along the Keldysh contour, J K dt — Ylri=± ^ = 
Sr;=± V I Expanding the exponentials in Ea. (|21ll we 
can write the second non-zero correction to S(t) as 




x ^(t)^"(0)5£(t 1 )«£(i 2 ))o. (22) 

The details of computing (S(t)) 1 are given in Appendix 
A. We obtain 



(Sit)), = -2e* 2 |r| 4 ( j£< 0<ta<t (/& - cosn(S, S„) + /^[costt^ + 5 a ) 2]} 

+.ftUfk cos tt((5 + S a ) + fig [cos 7r((5 +6 a )-l- cos 2tt^]} 

+ fir <0<t /«t[/« a «»' r (^ + ^)-/^cos7r(^-* tr )] 
+/§ [ft 2 cos 7r(<5 + 6„) - f 5 t g cos 7r(fy + 5 a )] 
+fto{fm cos7r((5 + 5 a ) + - cos^(<5 + 6„) - cos 2^]} 

+ CtT^ ~ f*o)[fk ™sn(5^ + S a ) ft costt^ - 8 r )] 

-ftUfk cos 7r(<5 + &„) + Fu\ [cos 2tt^ - 1 - cos tt(^ + 5 a )]}) (23) 

where 

& = i*nti - t 2 r\t - hi-^t - t 2 \-»\t 1 \-»\t 2 \->' cos[ Wo (i - h - 1 2 )} 

f? ti = \t - ^MltM* - t 2 \-^\h - t 2 \-^\h\-^ cos[wo(t + h- t 2 )] 

ft t2 = it-* 2 ntint-ti|- M ir' i iii-fei-it2|- M cos[wo(i+t2-ti)] ) (24) 
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fto t t = ftot t 2 ( u>0 - ¥ 0), and for the Pfaffian state 8 a — 5^ — 1/4. This reduces to: 
(S(t))t = -2e* 2 |r| 4 { r <0<t2< V/«t " fj)(fk ~ Vm) ~ fujw] 

ti<i 2 <0<t / -0<i 1 <t 2 <t 

\JtoJtt 2 JttiJto)' / \JtoJtt 2 itt\3tt-z t Ht%Ht\) ( \^°) 

l,t 2 Jti,t 2 

For comparison, we also provide the similar result for the case of a Laughlin state with filling factor obtained 
from the more general Pfaffian result by setting 6^ 2i>, and 5 a — > 0. 

<.ti<o<t 2 <t 

(S(t))i = -2e* 2 |rf{ / [2/ t 5 '(l - costtJO + 2/& (cos tt*, - 1) + // t ! 2 (2cos^ - 1 - cos2^)] 

J Ufa 

rtl < t 2 <0<t /.0<il<t 2 <t 

+ /5(l-coB27r*i)+ / // t ! 2 (l-cos27r^)}, (26) 




and 5; = 2*a 

A sketch of the integration domains is presented below: 




FIG. 3: Intervals of integration 



V. DISCUSSION AND RESULTS 

A few comments about (S(t))i are in order. Note that 
we focus only on t > 0, as (S(t)) is symmetric under 
t — ► — t. Also, when we study it numerically we need 
to subtract the corresponding constant (It) 2 , such that 
(S(t))i goes to zero for large time separations. The re- 
sult described above is for zero temperature, and can be 
translated to finite temperatures by a conformal trans- 
formation under which \t\ becomes sinh(|i|7rT)/7rT. 

We also note that we can expand cos[o;o(i+ei^i +£2^2)] 
as cos(ujot) cos[ujo( e iti + £2^2)] — sin(o;oi) sin[wo(ei^i + 
£2^2)]) where ei/2 = ±1- Thus, after performing the 
integrals over t\ and £2, (S(t))i can be written as 
cos(ui t)fi(.t) + sin(w i)/ 2 (i), where fi and / 2 will be os- 
cillatory decaying functions. The nice thing about this 
way of writing (S(t))i is that the contributions form the 



terms proportional to cos(woi), and sin(uJo)t in (S(t))i 
can be separated. In general the first type of terms will 
yield a Lorentzian peaked about loq when Fourier trans- 
formed, while the second type of terms will yield a deriva- 
tive of a Lorentzian type curve, also centered about the 
Josephson frequency ojq. 

Thus, we expect that, when Fourier transformed, 
each term in equation ll'MI) yields a superposition of 
Lorentzians and Lorentzian derivatives centered about 
the Josephson frequency ujq. Of course, the magnitude 
and the sign of each term depend on the interval of in- 
tegration, as well as on the corresponding statistical fac- 
tors. Also, the sharpness of each resulting Lorentzian 
peak/peak derivative will be determined by the temper- 
ature, as well as by the exponent of the corresponding 
power law. The exact shape of (S(lu))i will be thus de- 
termined by the interplay of many factors, and in general 
cannot be guessed without a detailed numerical evalua- 
tion of each term. 

Nevertheless, there a few observations one can make. 
The first is that in the case of a Laughlin state all the 
power laws involved in (<S(i))i come as ffo Ml>tt2 , while 

in the Pfaffian state some terms come as f^^* tt2 , while 

others come as combinations of f t Q tl t2 and //o tl t2 - A 
similar observation can also me made about the statisti- 
cal factors, in the Laughlin state all the statistical fac- 
tors involve Si , while in the Pfaffian state they will involve 
S^+Sa- and S^—Sa- . The sources of these difference are the 
non-abelian characteristics of the four point correlation 
function of the a operators. Specifically, as described in 
detail in Appendices A and B, the four point correlation 
function for the a operators has different forms, depend- 
ing on the order of the four times at which the correla- 
tion function is evaluated. Thus for each time ordering 
of t±,t2, and t, and thus for each domain of integration, 
each of the term to be integrated will have a different 
form. 

This will have a double effect, one on the relative 
coefficients, and one on the shape of each term. We 



can compare the relative magnitudes of each term. We 
note that for the Pfaffian case the integral is domi- 
nated by the terms involving f^* and /f t * on the interval 

< t\ < t- 2 < t. In particular, the term / t ** is responsible 
for the differences between the Laughlin and the Pfaffian 
cases and it vanishes in the Laughlin state, where the 
result is dominated by the term / tt * . 

The value of (S(u>))i, for a Pfaffian and Laughlin state 
(y = 1/5) are plotted below. We chose the values of 
the applied potential such that ljq = 1. Also we chose 
the temperature to be T = 0.1, and tunnelling strength 
|r| 2 = 0.5% (Pfaffian), and |r| 2 = 0.8% (Laughlin). The 
numerical integrals over t\ and t 2 are done using Math- 
ematica. The singularities at t = ii/2,0, and ti = t 2 
are avoided by stopping the integral at some S t = 0.001 
away from the singularities (this is equivalent to impos- 
ing a high energy hard cutoff). We could have imposed 
a soft cutoff as in the previous section by substituting 
it — ► e + it but the hard cutoff procedure is more trans- 
parent for the observation of statistics, and it is also eas- 
ier to implement numerically. Also, note the rounding 
of the Josephson singularities due to finite temperature 
effects. 
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FIG. 5: The total shot noise (S(u)) (in units of e* (It)) as a 
function of frequency (in units of u>o) in the Pfaffian state. 
The dotted line is the first order contribution to shot noise. 
The solid line is the shot noise coming from both the first and 
the second order terms. Note the shifting and the enhancing 
of the peak at the Josephson frequency as a result of including 
the higher order term, and the low frequency reduction in the 
shot noise. We set T = 0.1, and |r| 2 = 0.5% 



0.4 r 




FIG. 4: The second order correction to the frequency de- 
pendent shot noise (S(ui))i (in units of e* (It)), as a function 
of frequency (in units of uio) in the Pfaffian state. We set 
T = 0.1, and |r| 2 = 0.5% 

We note that the effect of the higher order correction 
in the noise is an enhancing and a shifting of the sin- 
gularity towards lower frequncics in the Pfaffian state 
and a flattening of the singularity in the Laughlin state. 
Also we note a decrease in the shot noise at zero and 
small frequencies, such that the ratio (S(u>))/(It) de- 
creases slightly in the Pfaffian state and in the Laughlin 
state with v = 1/3 and increases slightly for the Laugh- 
lin state with v — 1/5. This effect appears to be more 
significant in the Pfaffian state in the sense that a tun- 
neling strength L which produces comparable shifts of 
the Josephson frequency leads to a much larger suppres- 
sion of the zero-frequency shot noise in the case of the 
Pfaffian. 



FIG. 6: The second order correction to the frequency de- 
pendent shot noise (S(u))i (in units of el (It)), as a function 
of frequency (in units of uo) in the Laughlin state. We set 
T = 0.1, and |F| 2 = 0.8% 



This effect is not as dramatic as the underlying physics; 
once it is Fourier transformed and inserted into the ex- 
pression for the noise, the structure in the quasiparticle 
four-point function is masked. Hence, it might also be 
worth analyzing other setups, see for example the se- 
tups in Ref. ITllll2L A very similar calculation for abelian 
Laughlin and Jain states has been performed in a three 
terminal geometry^. The authors conclude that the 
dependence of shot noise in their setup also included 
Lorentzian peaks and peak derivatives, yielding contri- 
butions similar to our findings, and which could distin- 
guish between the Laughlin and Jain states. This setup 
is harder to achieve experimentally, but it has the ad- 
vantage that, while in our setup the relevant non-abelian 
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FIG. 7: The total shot noise (S(u)) (in units of e*(7t)) as a 
function of frequency (in units of wo) in the Laughlin state. 
The dotted line is the first order contribution to shot noise. 
The solid line is the shot noise coming from both the first and 
the second order terms. Note the shifting and the flattening 
of the peak at the Josephson frequency as a result of including 
the higher order term. We set T = 0.1, and jT| 2 = 0.8% 



physics comes in only as a second order correction, in a 
three terminal setup, it may be evident in the first non- 
zero order, which would make it easier to observe experi- 
mentally. In the next section, we make a few observations 
about the formalism needed for such a calculation. 

To sharpen the distinction between the Pfaffian state 
and other quantum Hall states (especially Abelian ones) , 
it would be useful to repeat our calculation for the Jain 
sequence, the (3,3,1) state, etc., and compare these re- 
sults with those given above. 



VI. MULTI-TERMINAL SETUP 

Since the effects of non- Abelian statistics first appear 
in fourth-order correlation functions, it would be advan- 
tageous to look at experimental setups in which the lead- 
ing processes occur at this order. Consider the setup de- 
picted in figure|Hl In order to describe tunneling between 
the different leads, we introduce multiple chiral Pfaffian 
edges <7j, ipi, 4>i, with i — 1,2, ... ,n. (These fields are 
all chiral, so we do not bother adding a subscript R.) 
Naively, the tunneling operator between edges i and j 
takes the form: 

= -IViOj-e'***-*^/ 2 + h.c. (27) 

However, this is not quite right because it treats the dif- 
ferent edges as completely independent. In reality, they 
are not quite independent because quasiparticles at the 
different edges must still have (non- Abelian!) braiding 
statistics. This can be accommodated by modifying 1|27|) 
to 

Hti = -TFjFj Ti t 3 a^e^-^/ 2 + h.c. (28) 



by introducing 'Klein factors' Fi and Tj, i — 1,2, ... ,n. 
These operators satisfy commutation relations which 
are 'reverse engineered' to give the correct quasiparti- 
cle statistics. The i^s account for the Abelian part and 
satisfy the relations^: 

F}Fi = 1 

FiFj = e^^FjFi (29) 

where pij — 0,±1 are chosen to ensure that tunneling 
operators commute if their tunneling paths do not inter- 
sect but acquire a phase ±iirg/4: upon commutation if 
the tunneling paths have intersection number ±1. (Re- 
member that g = 1/2 for a Pfaffian state at v = 5/2.) 

These relations were obtained by requiring that H^ n 
and i?tun have the same properties as anyonic quasipar- 
ticle trajectories (or Wilson lines). In the non- Abelian 
case, this generalizes to a skein relation. This skein rela- 
tion can only be satisfied with matrices, from which the 
non-Abelian nature of the statistics follows. Following 
the logic used in the ABelian case, but now generalized 
to this more complicated situation, we suggest that the 
correct statistics can be implemented by making the TjS 
satisfy the following relations. 

rf = 1 (30) 

To write down the commutation relations between differ- 
ent TjS, we group the leads into pairs. This can be done 
arbitrarily but, for convenience, we group them as (1, 2), 
(3, 4), . . ., (2j — 1, 2j), .... Each pair of leads has a gap- 
less fermionic level associated with it. The commutation 
relation between two r^s from the same pair is: 

T2j-lT2j = eT ff i T2jT2j—l (31) 

When a quasiparticle from one pair of leads hops onto 
a lead from a different pair, the two two-level systems 
are coupled. The commutation relation between two TjS 
from different pairs features a 4 x 4 matrix which acts on 
these two two-level systems 

^ = ^(_,^Tr nr k (32) 

The numbers pu are the same as for the Abelian Klein 
factors Fi. 

In the two-terminal case, which was our primary in- 
terest in this paper, these Klein factors should, strictly 
speaking be present, but they have no effect since the fac- 
tors associated with hopping from lead 1 to lead 2 and its 
conjugate commute, e.g. (tit 2 )(t 2 ti) = (t 2 ti)(tit 2 ) = 1. 
In the multi-lead case, Klein factors can have a dra- 
matic effect on correlation functions even in the Abelian 
case, where their commutation relations take a simpler 
f orm i2 J 22_ rp^g same j g j- rue j iere _ j n particular, we expect 

that some of the effects which we have found at fourth 
order in the two-terminal case will now appear at low- 
est order for current cross-correlations between different 
leads. This will be further investigated elsewhereiS. 
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I 

FIG. 8: A multi-lead setup. We treat the edge theory opera- 
tors at the different leads as independent, but must introduce 
Klein factors into our formalism in order to correctly account 
for quasiparticle statistics. 

VII. CONCLUSION 

We computed the shot noise in the tunnelling current 
between two edges of a quantum Hall fluid in a Pfaf- 
fian state. Specifically we analyzed the frequency depen- 
dence of the shot noise perturbatively in the tunnelling 
strength. We focused on the first two non-zero orders in 
perturbation theory. We found that, while the first non- 
zero order reveals some interesting physics, it does not 
encode any information about the non-abelian nature of 
the Pfaffian state. We also found that the second order 
in perturbation theory indeed contains information about 
the non-abelian physics. The effects of the non-abelian 



physics are manifest in the shot noise, in that the higher 
order correction has different features than, for exam- 
ple, a regular Laughlin state. Specifically, we found that 
the singularity in the shot noise at the Josephson fre- 
quency is shifted towards lower frequency and enhanced 
in the Pfaffian state and flattened in the Laughlin state. 
Also, the ratio of (S(t))/(I t ) at zero and small frequen- 
cies is decreased in the Pfaffian state and in the Laughlin 
state with v = 1/3, and increased in the Laughlin state 
with v = 1/5. Unfortunately, we cannot say if these fea- 
tures are entirely distinctive for a Pfaffian state and non- 
abelian statistics, and may appear also in other types of 
states, e.g. Jain, Halperin, 331. We hope however that 
our analysis and results will be a first step towards a thor- 
ough analysis of the effects of the non-abelian statistics 
in various other multi-terminal shot noise experiments, 
where more definite effects might be easier to isolate, 
and also will increase the interest to perform the shot 
noise experiments that would look for these sort of fea- 
tures. We must note that much of the relevant physics 
should come in at the Josephson frequency geV/2h in 
the range of GHz to THz, for applied voltages of order 
of microvolts to millivolts, which should be observable 
experimentally. 
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APPENDIX A: PERTURB ATIVE EXPANSION FOR THE CURRENT NOISE 

We will evaluate the second order correction to S(t): 

(Sit))! = ~ V(T K / dh [ rft 2 7"(i)/- r '(0)^(ti)5£(i 2 )}o. (Al) 
Z „=± J k Jk 

Here J K dt = X^=± V I^°oo ^ = Sr;=± V I dt. Also we have 

5C(t) = -rf(t)e- lbJot - F*f t e^ ot . (A2) 

and 

I t (t) = ie*TT{t)e- wat - ie*TT\t)e luat , (A3) 
where f(t) = a(t)e^ {t \ and e* = Thus we obtain for S(t): 

(S(t)) = -i( Je *) 2 |r| 4 EE E 61 f dt i I dt 2 e^ t+ ^-^+^ x 

V,Vi,V2 e=±lei=±le 2 =-l 

x(T K T e (t' ) )f" 1 (0- , ')f" 2 (^ 1 )f' e ( 1+ei+e2 )(^ 2 ))o. (A4) 
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We note that (S(t))i must be symmetric for t — > — t, thus we are only going to focus on t > 0. We also note that in 
this situation, various terms from the above expansion cancels for various ordering of the times t\, t 2 , and t. This is 
expected, due to causality, i.e., processes in which either one or both of two virtual quasiparticle tunnelling between 
edges happen after the measurement of (I t ) cannot contribute to fluctuations. In other words, only the situations 
for which at least one of the conditions t\ < < t and t 2 < < t holds will contribute. Noting the symmetry 
in exchanging t\ and t 2 , we conclude that we can write (S(t))i as a sum of three integrals over three domains of 
integration in the (ti,t 2 ) plane. The three domains of integration are depicted in Fig. |5| The intervals symmetric 
about the t\ = t 2 axis will provide an equal contribution, which is taken into account by an extra factor of 2. Given 




FIG. 9: Intervals of integration 



the fact that in each individual integration domain the ordering of the four times is well defined we can explicitly 
evaluate the Keldysh time-ordering operator. We thus obtain 

l-ti < o<t 2 <t 

(S{t)} = -e* 2 |r| 4 { / [F(0,t,t 2 ,h) + F(t,t2,0,ti) - F(0,t 2 ,t,ti) - F(t 2 ,t, 0,*i)] 

r ti < t 2 <o<t 



/ [F(0, t, t 2 , h) + F(t, 0, h, *i) - F{t 2 , 0, t, ti) - F{t 2 ,t, 0, h)} 

l-Q<tl < t 2 <t 

/ [F(0, t, t 2 , t x ) + F(t, t 2 ,t x , 0) - F(0, t 2 , t,h)- F(t 2 ,t, h,0)}} + h.c. (A5) 



where 



F(ai, o 2 , a 3 , a 4 ) = (a{a 1 )a(a 2 )a(a 3 )a{a4)){F ttl (ai, a 2 , a 3 , a 4 ) cos[uj (t + t x - t 2 )\ 
+F tt2 (ai, a 2 , 03, a 4 ) cos[u> (t + t 2 - h)] 

-F t0 (ai, o 2 , a 3 , a 4 ) cos[cj (i - * 2 - *i)]}. (A6) 
Here ai, a 2 , a 3 , a 4 are the possible permutations of of t\, t 2 , 0, t, and 

with t = 0, and s (0,t,ti,t 2 ) = (1,1,-1,-1), Sl (0,t,t l7 t 2 ) = (-1,1,1,-1), and s 2 (0, t, t u t 2 ) = (-1,1,-1,1). More 
explicitly F Ui is the four point correlation function of the four e*^ operators in which cf){t) and <f>(ti) appear with the 
same sign in the exponent. 

Using Wick's theorem we find that 

F tu (01,02,03,04) = |oi-o 2 |^ Sl(Ql)s * (Q2) |oi-03|^ Sl(Ql)Sl(Q3) |oi-04|^ Sl(Ql)Sl(Q4) x 

x \a 2 - a 3 \ s i' s ^ a2)si ^\a 2 - a 4 | <5 * Si(a2)s * (Q4) |o 3 - a4 \ s ^ {a3)St{ai) x 

x exp{i5^Tr[sgn(ai - o 2 )Sj(oi)sj(o 2 ) + sgn(a.\ - a 3 )s i {a 1 )s i {a 3 ) + sgn(ai - a 4 )s l (ai)s i (a 4 ) + 

+ sgn{a 2 - 03)51(02)5^03) + sgn(a 2 - 04)5,(02)8,(04) + sgn(a 3 - a4,)s l (a 3 )s i (a i )]/2}, (A8) 

where 8$ — g/2. 
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As detailed in Appendix B we also find the four point correlation function for the a operators starting from a CFT 
analysis of the correlation functions of the Ising model^I. 

(a{ ai )a{a 2 )a(a 3 )a{a 4 )) = F^(a x , a 2 , a 3 , a 4 )[0(1324) + 0(1342) + 0(2413) + 0(2431) + 0(3124) + 0(3142) 
+0(4213) + 0(4231)] 

+F? 3 (ai, a 2 , as, a 4 ) [0(1234) + 0(1432) + 0(2143) + 0(2341) + 0(3214) + 0(3412) + 0(4123) + 0(4321)] 
+F[ 4 (ax,a 2 ,a 3 , a 4 )[0(1243) + 0(1423) + 0(2134) + 0(2314) + 0(3241) + 0(3421) + 0(4132) + 0(4312)], (A9) 

where 6{abcd) = 1 for a a > otb > a c > ay and is zero otherwise. Also 



Ff 2 (ai, a 2 , a 3 , a 4 ) = \ot\ - a 2 \ 5 " \a 3 - a 4 | 5 " \a\ - a 3 \ 5<r |ai-a 4 | s " \a 2 - a 3 \ 5 "|a 2 -a 4 ' 



x 



x e i^S„[sgn(ax—a2)+sgn{aa—a4)—agn(ax—a 3 )—sgn(ax—a4,)—sgn(a2—a3)—sgn{a.2—a4,)]/2 (AlO) 

and similarly 



Ff 3 (a 1 ,a 2 ,a 3 ,a 4 ) = \a x - a 3 \ S " \a 2 ~ a 4 \ S " \a x - a 2 \ s<r |ai-a 4 | S "\a 2 -a 3 \ s "\a 3 -a 4 \ 



x 



x e iirS a [sgn(cti— a 3 )+sgn{ai2 — 0:4)— sgn(ax— 012)— sgn(eti— 014)— sff«(«2— "3)— sgn(a 3 — 014,)]/% (All) 



and 



F[ i (ai,a 2 ,a 3 ,a 4 ,) = \a x - a 4 \ s ° \a 2 - a 3 \ S " \a x - a 3 \ 5 "\ai-a 2 \ S "\a 3 -a A \ 5 °\a 2 -a 4 \ S " X 

x e inS a -[sgn{ai-at)+sgn(a2—a3)-sgn(ai-a 3 )-sgn{ai-a2)—sgn(a3—ai4)-sgn(a2-cit)]/2 (A12) 

with 8 a = 1/4. 

After some tedious but straightforward algebra one gets: 

(S(t)>! = -2 e * 2 |lf( j£< 0<t><t (/& - coBTrfo - ^) + //<f [costt^ + «5 CT ) - 2]} 

+ft 2 {ftt2 cos tt(^ + 6„) + /& [cos tt(^ +5 a )-l- cos 2tt^]} 

+ftt 2 [fk 2 costt(^ + <5 ff ) - //(f costt((5 + <5 CT )] 
+/»{/« cos7r((5 + <5 CT ) + - cos^(<5 + <5 CT ) - cos 27^]} 
+ C t ' 2 1<t2< * (/ t 5 / 1 -/^)[^cos7r(^ + ^)-/^co S7 r(^-5 (T )] 

{/& cos tt(^ + S a ) + /& [cos 2tt^ - 1 - cos tt(^ + <5 CT )] }) (A13) 

where 

/& = |*n*i - i 2 nt - tij-^l* - tsr^ltil-^ltai-^ cos[o;oC* - *i - * 2 )] 
= |t-tini 2 ntMt-t 2 |-' i |^-t2Mti|- M cos[a;o(t + ii-t2)] 

= |t-t a HtiHt-*i|"' , |t|"' , |*i-ta|" /i |t2|"' , co8[wo(* + *2-ti)], (AM) 

fto,t u t 2 = /totit 2 ( a; o -* 0), and for the Pfaffian state £ CT = <^ = 1/4. 

For comparison, we also provide the similar result for the case of a Laughlin state with filling factor obtained 
from the above equation by setting 8$ — > 2v and 8 a — > 0: 

rtl < 0<t 2 <t 

(S(t)) 1 = -2e* 2 |lf{ / [2/*j(l - costtJO + 2/& (cos ttJ, - 1) + (2 cos ttJ, - 1 - cos2^)] 

Jtl,t2 

tl < t 2 <0<t pQ<tl < t 2 <t 

/5(l-coB27r*,)+ / /// 2 (l-cos2^)}, (A15) 

1 1 ,*2 ** tl 5^2 

where 5/ = 2is. 
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APPENDIX B: FOUR POINT CORRELATION FUNCTION FOR THE a OPERATORS 



In order to derive the real time four point correlation function for the Ising a operators, we start from the imaginary 
time conformal field theory results^!: 



(a(zi)a(z 2 )a(z 3 )a{z4,)) = 



Z13Z24 



^14^23^12^34 



1/4 h+\x\+\i-X\ 



(Bl) 



where z — t + ix, and Zij — Zi — zj. We chose to denote X — z\ 2 z 3 ^j z\ 3 z 2 ^. Since all our correlations will involve 
operators at the tunnelling point we will set all x^s to zero. In this limit we can evaluate the four point correlation 
functions and we note that 



(a(z 1 )a(z 2 )a(z 3 )a(z i )) = 



Z13Z24 



Z14Z23Z12Z34 
Z13Z24 



1/4 



for < X < 1, 



Z14Z23Z12Z34 
Z13Z24 



1/4 



1/4 



VT~x = 



Z14Z23 



Z13Z24Z12Z34 



1/4 



Z14Z23Z12Z34 

By analytical continuation tj — > iai to real time this yields: 



Vx = 



Z12Z34 



^13^24223^14 



1/4 



, for X < 0, 



for X > 1, 



(B2) 



(a(ax)a(a 2 )a(a 3 )a(a 4 )) = F? 2 (a u a 2 ,a 3 , a 4 ) [0(1324) + 0(1342) + 0(2413) + 0(2431) + 0(3124) + 0(3142) 
+0(4213) + 0(4231)] 

+F[ 3 (a 1 ,a 2 ,a 3 , a 4 )[0(1234) + 0(1432) + 0(2143) + 0(2341) + 0(3214) + 0(3412) + 0(4123) + 0(4321)] 
+F[ i (ai,a 2 ,a 3 , a 4 ) [0(1243) + 0(1423) + 0(2134) + 0(2314) + 0(3241) + 0(3421) + 0(4132) + 0(4312)], (B3) 



where 6{abcd) = 1 for a a > at, > a c > ay and is zero otherwise. Also 
Ff 2 (ai,a 2 ,a 3 , ar 4 ) 



x e 



on - Oi 2 \ s '\a 3 - a A \ 6 "\ ai - a 3 | 5<T K - a 4 | A<T |a2-a 3 | s "\a 2 -a A \ 6 " x 

iTr&a [sgn(pL\ — OL2)~l-sgn(oL3— ol^) — sgn^OL\ — as) — sgn(ai — a&) — sgn(a2 — as) — sgn{a2 — a&)]/2 



and similarly 



and 



F[ 3 (ai,a 2 ,a 3 , ar 4 ) 



F[ 4 (ai,a 2 ,a 3 , ar 4 ) 



x e 



\oti — as\ s "\a2 — a 4 | 5 "|ai — a 2 \ ^lai — a 4 | 5 " \a 2 — a 3 \ s "\a 3 — a 4 | s " x 

iTT&a [sgn{ai — as) J rsgn(a2—a&) — sgn(ai — CK2) — sgn(a± — a&) — sgn(a2 — as) — sgn(as — 04)]/ 2 



x e 



\ot\ - ot4\ 5 °\a2 ~ a 3 \ s °\ai - a 3 \ 5 "\a 1 -a 2 \ ^ |a 3 - a 4 | ^ |a 2 - a 4 | S " X 

i7r5^[sgn(ai— a 4 )+sgn(a2— 0:3)— S9«(ai— "3)— s 9»(«i— "2)— sff«(a3— «4)— S9«(a2— a*)]/2 



(B4) 



(B5) 



(B6) 



with 5o- = 1/4. 

We can try to understand the physics of this result; one notes that F° 2 ^ 13 ^ 14 are the three cross ratios appearing 
in the fourth order correlation functions. In a fourth order correlation function for an abelian state, they usually 
appear in quasi-symmetric combinations, such that by exchanging two of the times the most significant effect is the 
apparition of various phase factors, characteristic of fractional statistics. In the non-abelian case however, we see that 
the physics is very different, in that when we exchange two of the times, not only one acquires phase factors, but one 
can also change the form of the correlation function from one of -F 1 ^ / 13 / 14 cross ratios to another. 
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